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Abstract 

We study chiral symmetry breaking in QCD-like gauge theories introduc- 
ing a confining effective propagator, as proposed recently by Cornwall, and 
considering the effect of dynamical gauge boson mass generation. The ef- 
fective confining propagator has the form l/(/c 2 + m 2 ) 2 and we study the 
bifurcation equation finding limits on the parameter m below which a satis- 
factory fermion mass solution is generated. Considering the evidences that 
the coupling constant and the gauge boson propagator are damped in the in- 
frared, due to the presence of dynamically massive gauge bosons, the major 
part of the chiral breaking is mostly due to the confining propagator. We 
study the asymptotic behavior of the gap equation containing confinement 
and massive gauge boson exchange, and find that the symmetry breaking 
can be approximated at some extent by an effective four-fermion interaction 
generated by the confining propagator. We compute some QCD chiral pa- 
rameters as a function of m, finding values compatible with the experimental 
data. Within this approach we expect that lattice simulations should not 
see large differences between the confinement and chiral symmetry break- 
ing scales independent of the fermionic representation and we find a simple 
approximate relation between the fermion condensate and dynamical mass 
for a given representation as a function of the parameters appearing in the 
effective confining propagator. 
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1. Introduction 

QCD has two main properties: the chiral symmetry breaking (CSB) and 
confinement of quarks and gluons. Both phenomena are related to the non- 
perturbative infrared (IR) dynamics. The dynamical generation of quarks 
masses leading to CSB can be observed through the study of QCD Green 
functions, what has been extensively performed in the last years by means of 
Schwinger- Dyson equations (SDE) [1]. Green functions can also be studied 
through simulations of gauge theories on the lattice, and, in particular, a 
great improvement in the understanding of dynamical gauge boson mass 
generation in these theories was obtained recently as pointed out in Ref . jij • In 
what concerns gauge boson mass generation there is a nice agreement between 
the lattice computations and the solutions of Schwinger- Dyson equations 
(3], Q . The consistency between these two different approaches, strengthened 
by the phenomenological consequences (if, reinforces the robust picture of 
auge boson mass generation formulated by Cornwall many years ago 0, 0, 



8|, |9|, [10|, 111] . In this scenario the dynamical gluon mass induce vortices in 
the theory and these may be responsible for confinement. 

The study of CSB through lattice simulations and the SDE approach 
have not come to the same firm ground as in the case of gauge boson mass 
generation. First, early SDE studies of CSB have shown a minute or absence 
of quark mass generation in the presence of dynamical gluon mass (m g (k 2 )) 



12l Il3l . Il4fl . Secondly, in lattice simulations it seemed that CSB and con- 



finement were triggered by the same mechanism (at least for SU(2) theories 



15|,|16(). A recent result about CSB in SU(3) gauge theory indicates that we 
still need more simulations to confirm whether or not confinement and CSB 
are intimately connected, but there is a clear sign of a connection between 



CSB and the string tension [17|. From the point of view of gap equations 
we may say that something is missing in the solution of the fermionic SDE, 
and the most plausible possibility is confinement at least for fermions in the 
fundamental representation, although it may not be essential for CSB when 
fermions are in higher dimensional representations (e.g. "quarks" in the ad- 
joint representation) [l8[. Two recent papers tackle this problem in different 
ways: one introduces a confining effective propagator into the gap equation 
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19[ and the other introduces a modified quark-gluon vertex containing infor- 



mation about the ghost sector as well as makes use of lattice propagators [20 
Both studies are able to explain CSB for fundamental and adjoint quarks. 
One of them indicate that confinement may be an essential ingredient in the 



gap equation [191 ] . The other contains a set of effects which include the use 
of a gluon propagator obtained from the lattice simulations and a complete 
vertex function that is enhanced in the infrared 2^ . 



If CSB is a phenomena linked to confinement, as the lattice results seem 
to indicate, we may ask how the confinement information is embodied into 
the lattice gluon propagator used in the fermionic gap equation in order to 
generate quark masses. Anyhow, some gluon confinement effect ought to 
be present in the lattice propagator. However it does not show the linear 
confining potential felt by quarks that has also been observed in the lattice 
simulations [21] . The linear potential that we are referring to is the successful 
phenomenological quark model potential given by 

V F (r) = K F r-^, (1) 
3 r 

where the confining first term is linear with the distance and proportional 
to the string tension Kp. The second term, that is of order a s , the strong 
coupling constant, describes the one gluon exchange contribution. The fact 
that the lattice gluon propagator does not reproduce the confining part of 



this potential has been observed recently in Ref.[22|. The potential between 
static quark charges is related to the Fourier transform of the time-time 
component of the full gluon propagator in the following way 

= J dV s (q 2 )A 00 (q) exp^ r , (2) 

where Gi is the Casimir eigenvalue of the fundamental representation, the 
bold terms, q and r, are 3- vectors. A o(q) is the zero-zero component of the 
gluon propagator in the momentum configuration. As noticed in Ref. 22j the 
linear confining term of the potential (Kpr) cannot be obtained from the 
gluon propagator determined in the lattice or from the gluonic SDE. Only 
the Fourier transform of a l/q A type of the product coupling®propagator will 
generate a linear term in the quark potential, although no fundamental field 
presents such behavior. 



The proposal of Ref.[l9j is attractive because it includes the necessary 



confining force into the gap equation, maintaining covariance and solves the 
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problem discussed in the previous paragraph. It is a simple procedure based 
on critical properties of confinement and, as we will discuss, allows for a 
phenomenologically consistent discussion of CSB in gauge theories, providing 
a model that may be tested through lattice calculations. The main idea is 
that entropic effects, that can be associated to an area law for the Wilson 
quark loop, can be taken into account by a confining effective propagator 
of the form l/(p 2 + m 2 ) 2 with a finite parameter m oc Kp/M(0), where 
M(0) is the zero-momentum value of the running quark mass and Kp is 
the string tension. Using this confining effective propagator, which is to be 
understood as the result of a collective effect not related to the propagator 
of a fundamental gluon field, and responsible for quark confinement. Given 
this proposal within its own reasons, we aim to investigate its consequences 
for CSB phenomenology in the QCD case. The new gap equation is able to 
reproduce all QCD chiral symmetry brea king parameters. 



We may say that the papers of Ref. [20] and 19| have complementary 



ideas and may even indicate a possible mixed picture of the CSB mechanism. 



In Ref. 20] it was noticed that even a quite complete treatment of the vertex 



function and gauge dependence of the quark SDE does not allow for CSB in 
agreement with experiment. A gluon propagator determined in the lattice, 
that is a little bit less damped in the intermediate infrared region when 
compared to the one obtained with the gluon SDE, is necessary, as well as 
some strong ghost effects into the quark-gluon vertex. However, it is not 
clear that introducing some complicated ghost physics into the one-gluon 
gap equation would lead to confining effects typical of center vortices. It 
is also questionable that the gluon propagators found on the lattice will be 
able to evidence confinement in, for example, a Wilson loop. If so much 
confinement is present in the one-gluon-exchange approach it would not be 
simple to explain how the pion ends up massless. Finally, it is not through 
open quark lines, as all gap equations use, that we shall observe confining 
effects, but through closed loops. Hence there is certain appeal on the explicit 
introduction of confining objects, such as center vortices. This is exactly what 



Cornwall's proposal does [19(, introducing an effective confining propagator 
whose effect is to add an extra strength to the massive gluon exchange in 
the infrared, with the advantage of having a simple formula to model the IR 
QCD behavior. We also recall that it has been argued that the one-gluon 
massive exchange and the confining gluon propagators act differently when 
quarks are in the fundamental and adjoint representations 18| . therefore we 



have a precise way to investigate the possible differences in the CSB and 
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confinement transitions for fermions in different representations. We will 
argue that CSB is a direct consequence of confinement in this specific model. 
The article is organized as follows. In Sect. II, we describe briefly the 



model of Ref.[19|, i.e. how confinement is introduced in the gap equation. 
In Sect. Ill, we study the bifurcation equation of the full gap equation. It 
is verified that we have maximum m values in order to generate reasonable 
values for the dynamical quark mass. In Sect. IV, we present an analysis of 
the asymptotic behavior of the gap equation. The idea was to investigate how 
confinement could affect the asymptotic behavior of the fermionic self-energy. 
It must be said that the CSB mechanism is not only important for QCD, but 
for any non-Abelian gauge theory. In particular, it has consequences for 
the Standard Model dynamical gauge symmetry breaking, or the so called 
Technicolor models. We find that if the confinement effective propagator were 
restricted to a limited region of momenta the self-energy is well approximated 
by an effective four-fermion interaction. In Sect. V we compute quantities 
like the pion decay constant and the quark condensate verifying that the 
model predicts reasonable values for the QCD chiral parameters. In Sect. VI 
we determine an approximate relation between the fermion condensate and 
the dynamical mass for a given representation as a function of the parameters 
appearing in the confining effective propagator. These quantities can be 
studied through lattice simulations and may provide a test for this CSB 
mechanism. In Sect. VII, we draw our conclusions. 



2. Introducing confinement into the gap equation 

The possible relation between confinement and chiral symmetry breaking 
for quarks in the fundamental representation is an issue discussed several 
times in the literature in many different contexts as pointed out in Re fjl8 ]. 
With the evidence that gluons acquire a dynamically generated mass @, Uli|, 
i.e. a momentum dependent mass m g (k 2 ), it also becomes clear that the 
standard quark gap equation with dynamically massive gluons does not have 
enough strength to generate quark masses. Therefore, if CSB and confine- 
ment are triggered by the same mechanism, the main quest in this subject is 
how to introduce confinement into the quark SDE. 



The proposal of Ref.[19| is that confinement should be introduced into 



the gap equation through the following effective propagator, meaning that it 
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is not at all related to the propagation of a standard quantum field: 

8irK F 
[k 2 + m 2 ) 2 



D% f (k) = &»D eff (k); D eff (k) = j^^A^ , (3) 



In the m — V limit we obtain the standard effective propagator %irKp5^ u /k A , 
that yields approximately an area law for the Wilson loop. We must neces- 
sarily have a finite m^fl value due to entropic reasons as demonstrated in 



Ref . [19[ , and it is related to the dynamical quark mass, as required by gauge 
invariance. Moreover, the Abelian gauge invariance of this effective propa- 
gator must appear in the quark action obtained by integrating over quark 



world lines that will imply a area- law action [19(. With the inverse fermionic 
propagator written as 

S-\p) = p'A(p 2 ) + B(p 2 ) , 

using the approximation A(p 2 ) = 1, which implies that B(p 2 ) turns out to 
be identical to the dynamical mass M(p 2 ), we can see that the effective 
propagator of Eq. (j3J) leads to an Abelian gluon gap equation equal to 

M ^=(kJ d ' kD <>>^- k) ^m' (4) 

where M c (p 2 ) is the running dynamical quark mass generated by the confining 
propagator, so we end up with 

M h 2 ) ~ 1 fd'k 8nKF 4M(fc2) (5) 

c[P ] (2vr) 4 J [(p - k) 2 + m 2 } 2 k 2 + M 2 (k 2 ) ' 1 } 

Note that there is an interplay between the parameters Kp, m and m„ as 



extensively discussed in Ref.|19|. and we will comment about their values, 
but it must be said that they should be related among themselves because 
in the end all of them will be proportional to the QCD mass scale Aqcd- 

Eq.(jSJ) is not the whole story, since it was also recognized that massive 
one-dressed-gluon exchange may induce a quite small CSB for fermions in 
the fundamental representation as well as a larger symmetry breaking for 



fermions in higher dimensional representations [18|, [20J, which, in the Landau 
gauge, is given by 



( 2) = [ d , g 2 (p-k)3M(k 2 ) 

W J (2vr) 4 J [{p-k) 2 + m 2 {p-k)][k 2 + M 2 {k 2 )} 
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where C2 is the quark Casimir eigenvalue and g 2 is the effective charge 

1 



f(k 2 ) 



bH(k 2 + 4m 2 g )/A 2 QCD ] 



(7) 



where b = (UN —2n f) / 48tt 2 for the SU(N) group with n/ flavors. For quarks 
in the fundamental representation C2 = 4/3 and with the phenomeno logically 



preferred value m g ps 2Aq C d ~ 500 — 600 MeV |5[|, we see that this charge's 



value at the infrared fixed point |25|, |26| (a s (Q) = 5 r2 (0)/47r) is of order 0.5, 
as shown in Fig.flT]), while it should be at least a factor 2 larger to trigger 
CSB. 



Running coupling 

m =500 MeV and n =3 



1E-3 0,01 0,1 1 10 100 1000 10000 



Figure 1: The running coupling constant for QCD with dynamically massive gluons. 



For "quarks" in the adjoint representation the Casimir eigenvalue is approx- 
imately a factor two larger and compensates the small effective coupling 



18l . |20| . The ultimate gap equation may be quite sophisticated but it can be 



well modeled by the sum of the "confining" plus massive one-gluon exchange 



19|, which, in the Abelian gluon approximation, is given by 



M(p 2 



4M{k 2 



(2k) 



d A k 



k 2 + M 2 (k 2 ) 
g 2 (p- k)3M{k 2 ) 



[(p - k) 2 + m 2 Jp- k)} [k 2 + M 2 (k 2 )} ' 



(8) 
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where M(p 2 ) = M c (p 2 ) + Mi g (p 2 ) is the dynamical quark mass generated by 
the confining and one-dressed-gluon propagators. This last equation is the 
basic one that we shall explore in this work. We could say that this equa- 
tion resembles, in a different context, what we have in the phenomeno logical 
quark model potential described in Eq.flT]), i.e. a part that is responsible by 
confinement (generating an approximately linear term proportional to the 
string tension), and the second term, that is of order a s , describing the one- 
gluon exchange contribution. Thereof, confining part of Eq. (JS]) is a reasonable 
phenomenological way to study CSB taking into account the effective area 
law. 

The solution of the confining gap equation (Eq.flS])) was discussed by 



Cornwall [19[, and it was observed that the confining propagator gener- 
ates CSB represented by a dynamical quark mass that is of order M ~ 
m ~ y 1 Kpjix w Aqcd- This dynamical quark mass has a fast large mo- 
mentum falloff (oc 1/p 4 ), which results from a gap equation basically dom- 
inated by the infrared region of the propagator. On the other hand the 
1-gluon gap equation was extensively studied in the literature and it is 
known that it generates CSB only above a certain critical coupling equal 
to a c (0) = (g 2 /4n) > (7r/3C 2 ), which, as discussed in the previous para- 
graph, is not achieved by the effective coupling of Eq.Q. Notice that the 
suggestion of a hypothetical Casimir scaling law for CSB with fermions in 
the representation R 23|, \2A 



a s C 2 (R) « 0(1) , (9) 

appears as a consequence of the SDE equation without the existence of a 
dynamical gauge boson mass. 

It is easy to verify that Eq.flSJ) has nontrivial solutions. Its critical behav- 
ior can be inferred from the derivative of Eq.((6]) with respect to m 2 , where 
we can replace ?>g 2 Ci by %i\Kp (in the constant coupling limit), linearize 
with the substitution of M 2 (k 2 ) in the denominator by M 2 , and evaluate at 



m g = m: 



oc M cOT . (10) 



dm 2 



Eq.(|I]) has already been solved with different approximations |l2l. Il3l . Il9 
and has solutions only for large g 2 C*2 values as noticed in Ref. [19] . The 
condition for Eq. fflQ]) to have a nontrivial solution (and consequently Eq.flS}) 
is the same one that Eq.([n]) has to obey. Considering the damped confining 
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propagator that leads to an effective coupling 87i(K F /rn 2 ) (that is now in the 
place of 3g 2 C 2 ), and the fact that Kp ^ (3 — 4)m 2 (l9| . we verify, just by 
transposing the results of Ref. 12, 



19l 1271 . |28| (see, particularly, the simple 



analysis shown in section 3 of Ref. 13] ), that chiral symmetry will always be 



broken with the confining propagator. Therefore, the above procedure leads 
to the solutions of Eq. (jSJ) that can be represented by a linear combination of 
the following expressions: 

7±+l 



Mf{p) « M6{m 2 -p l ) + M 



2m 2 



p2 _|_ m 2 



#(p — m 



obtained from the derivative of the solutions of Eq.®, where 9 is the step 
function and 



7± = \ [1 =F (1 



4A) 1/2 ] , \ = K F /2ix 



The value of M is determined by the boundary condition at p 2 = 



M = M c (p 2 = 0) 



2K 



71 



dk 2 



k 2 M{k 2 



(k 2 + m 2 ) 2 (k 2 + M 2 ) ' 



(12) 



It seems that the confining effective propagator for the Kp and m values 
discussed above is able to describe the chiral s ym metry breaking. However 
the full CSB problem, as discussed by Cornwall [l9| , includes the effect of the 
dressed 1-gluon exchange, and the critical behavior for the onset of dynamical 
quark masses in this case will be discussed in the next section. It is possible 
that only the sum of these effects may explain the lattice results for CSB 



17 . 



3. Critical behavior of the complete gap equation 



We can discover some aspects about Eq.([5I) critical behavior examining 
its bifurcation equation. To verify at what point the nontrivial solution of 
Eq.® bifurcates away from trivial solution, it is sufficient to consider the 



linearized version of that equation [27j, |28[| . We will deviate from the standard 
bifurcation theory proceeding as in Ref. 29|, and instead of substituting k 2 + 
M 2 (k 2 ) by k 2 in the denominators of Eq. (15|) , we will replace this term by 
k 2 + 5M 2 (0) and define the IR value of the dynamical quark mass (M) by 
the normalization condition 



SM(0) = M 
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following these steps we arrive at our bifurcation equation 



(2tt) 4 J [(p - k) 2 + m 2 } 2 k 2 + M 2 

^ f t- k)3 ^P „■ (13) 



(2tt) 4 7 [(p - k) 2 + m 2 ,] [P + m 2 } 

Our main intention in this section is to verify the gross critical behavior of 
the gap equation with these infrared finite propagators, this is the reason to 
select a bare massive gluon in Eq. (IT5|) . This equation is a standard Fredholm 
equation with a positive kernel, and, requiring 5M(p 2 ) to belong to L 2 , the 
spectrum is discrete with a smallest value for the "effective coupling" Kp/m 2 
and the 1-gluon exchange coupling g 2 /4:iT such that we have the trivial solu- 
tion SM(p 2 ) = for values of these couplings smaller than a certain critical 
value, and the nontrivial one if their values are larger than this same critical 
value. 

We can still make some simplifier approximations before estimating the 
critical behavior of Eq. (I13j) . making the following substitutions 

« [iP - k) 2 /A 2 QCD ] = -+ 9(p 2 -k 2 )a(p 2 )+9(k 2 -p 2 )a(k 2 ) , (14) 

and 

[p — k) A + rrig p A + k z + 

which is known as the angle approximation, and introduces an error of about 
10% in the calculation [30|. If we also define the variables x = p 2 /M 2 , y = 
k 2 /M 2 , k = m 2 /M 2 } e = m 2 g /M 2 } p = A 2 QCD /M 2 and f(p 2 ) = 5M(p 2 )/M } 
we obtain 



1 



7T 



A 2 /M 2 



f(x) = - dyK(x,y)f(y) , (16) 



o 



where we introduced an ultraviolet cutoff (A) and the kernel K(x, y) is equal 
to 



K(x,y) 



y 



(y + i) 



M 2 {y + n) 2 16vr (y + e) ' {U ' 



+ (™L 1 )9(x-y) 
' M 2 (x + k) 2 + 16tt (x + 1 1 y) 



(17) 
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The kernel K is square integrable 



1 11 ' ' ' V (y + 1) 2 \M 2 (x + k) 2 16vr (x + e) 



o Jo 



f IM2 . f IM2 , V 2 f2K f 1 , 3C 2 ^ 2 (y) \ 2 



J x 



therefore Eq.( fl6|) has a nontrivial L 2 solution for its coupling on a point set. 
The smallest eigenvalue (which, as we shall see, is related to the value of 
Kp/M 2 ) for which Eq. (fl6|) has a nontrivial square integrable solution, is the 
first bifurcation of the nonlinear equation, and satisfies 

-l|iq = i. (19) 

71 

Note that the kernel K contains the sum of two contributions, that we may 
denote K c and Ki g , corresponding to the confining and dressed 1-gluon ex- 
change propagator. From the triangle inequality we know that 

\\K C + K lg \\ < \\K C \\ + \\K lg \\ , (20) 

from which we can recover the early results for the 1-gluon exchange, i.e. if 
the gluon mass is not introduced into the gap equation and the confining 
propagator is neglected, the limit of Eq. (JT9~j) is obtained only due to the 
Ki g contribution with a critical coupling constant a c ~ 0(1), but the K c 
contribution is necessary when considering dynamically massive gluons and 
quarks in the fundamental representation. 



The bifurcation condition described by Eq. (fl9|) is depicted in Fig. ([2]) in 
the case of n f = 3, Kqcd = 300MeV and K F = 0.18GeV 2 , where it is shown 
a dashed (blue) curve obtained with m g = 600MeV, a dot-dashed (red) curve 
with trig = 650MeV and a solid (black) curve with m g = 700MeV. Each point 
of these curves indicate the bifurcation point for a given m value generating 
a dynamical quark mass M. It should be noticed that there is a maximum 
m value above which there is no CSB, and this maximum value does not 
vary much as we change the values of the dynamical gluon mass. It is also 
interesting to verify that CSB also receives contributions from the massive 
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O 0.25 - 



0.10 0.15 
M(GcV) 



Figure 2: The criticality condition given by Eq. ([TO|) for the kernel of Eq. ([T7| is plotted 
in the case of n f = 3, A QC d = 300MeV and K F = 0.18GeV 2 . The dashed (blue) curve 
was obtained with m g — 600MeV, the dot-dashed (red) curve with m g = 650MeV and the 
solid (black) curve with m g — 700MeV. 



gluon term, and this is the reason for the differences between the curves, 
though the massive gluon (or 1-gluon exchange) generates a minute mass, 
as already observed many years ago in Ref . 12|, ll3(] . As the dynamical gluon 
mass is decreased we can observe a small increase in the maximum value 
of the m parameter, due to the 1-gluon exchange increasing contribution to 
the gap equation. When the bifurcation condition is computed it is possible 
to verify a larger instability in the numerical procedure as we go to larger 
values of m and m g (for large M values), which is due to the fact that we 
have two contributions to the chiral breaking, but with one of the kernels 
contributing much more to the breaking than the other. This fact is obvious 
if we compare the different propagators that we are dealing with. These 
propagators are shown in Fig. ([3]) and Fig(HJ) for typical mass parameters 
that we are discussing in this work. We can see a huge difference between 
the two propagators (O(10 3 )). Of course, the confining effective propagator 
should not be regarded as the actual gluon propagator, but must be seen as 
a collective effect that produces confinement^ with an area law for the quark 
action and appropriate entropic properties [l9(. The difference between these 
contributions are responsible for the delicate convergence of the bifurcation 
condition. 
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It must be also noticed that the confining effective propagator has most of 
its effect concentrated in a momentum region below (9(100) MeV. Actually, 
as we shall discuss ahead, most of the chiral symmetry breaking will occur 
due to the low momentum region of the gap equation, what is consistent 
with lattice observations that the relevant momentum component of gluons 
for CSB is exactly this region 3l|, |32|, |33|. This fact may be contrasted 



with the results of Ref.[20], where a larger part of the CSB comes from an 
intermediate region of 0(1) GeV. 




1E-3 0,01 0,1 1 10 100 1000 



k 2 [GeV 2 ] 



Figure 3: The massive gluon propagator with n/ = 3 and Aqcd = 300 MeV. 



4. The asymptotic behavior of the gap equation 

The asymptotic behavior of the complete gap equation [Eq.fjH])] can be 
obtained from the linearized Eq. (I16[) with the kernel given by Eq.( lT7|) in 
a procedure identical to the one performed by Takeuchi [34| and Kondo, 



Shuto and Yamawaki [35] . In these references the QCD CSB problem was 
solved considering a quark Schwinger- Dyson equation also with a two kernel 
contribution: One due to an effective four-fermion interaction and another 
one due to a perturbative gluon exchange. 

The integral equation (TIB]) with the kernel described in Eq. (ll7j) can be 
transformed into a differential equation and solved with appropriate bound- 
ary conditions. In order to still simplify the calculation we assume 

K = e = p=l, (21) 
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Figure 4: The confining effective propagator of Eq.([3]) with Kp = 0.18 GeV 2 and different 
m values. Note the 0(1O 3 ) strength difference when compared to Fig. ((3]). 



and neglect the running of the coupling constant in Eq. fj!7j) . The fact that 
we consider all masses equal to M, as implied by Eq. (j2ip . shall not modify 
the asymptotic behavior of the quark self-energy, and the assumption of a 
constant coupling will introduce sub-leading or logarithmic corrections to 
the asymptotic solution. After these approximations we obtain the following 
differential equation 

(x + l) 3 [2a 1 (x+ 1) + a 2 (x + l) 2 ]f"(x) + (x + lf[6a 1 + 2a 2 (x + l)]f(x) 
+ x[2at + a 2 (x + l)] 2 f(x) = , (22) 



where 



2K f 3C 2 g 
ai = ^P ' a2 = T6V ' (23) 
with the infrared and ultraviolet boundary conditions given respectively by 

/(*)U> = 1 » ( 24 ) 

The asymptotic solution of the linear second-order differential equation, 
with a singularity at infinity, is obtained as a linear combination of two 
independent solutions, f(x) = b\ f+(x) + b2f-(x), which can be obtained by 
applying the expansion method [36|. Eq. (l22]) can be put in the form f"(x) + 
p(x)f + q(x)f(x) = 0, where p(x) and q(x) can be expanded in convergent 



14 



series of the form p(x) = (1/x) J2^=oPs x s an d q( x ) — (V^ 2 ) Yl^Lol^ S ^ or 
large x. The asymptotic solutions are of the form 



= ^£§> ( 25 ) 

s=0 

where a, are the roots of 

oci(oii + 1) - Poca + qo = . (26) 

With the leading coefficients of the differential equation given by po = 2 and 
go = 0,2 we obtain two roots 



«± = 2 = — • (2?) 

Defining a new variable 

w±l , , 

7± = — , (28) 

we can write the asymptotic solution as 

OO _|_ oo _ 

/(*) = hx^ + W~ £ | ( 29 ) 

The substitution of Eq.()25l). with 7± given by Eq. (|28l) . into the differential 
equation (1221 allows the determination of the c\ coefficients and the recursion 
formula satisfied by them. 

Considering only the "+" solution displayed in Eq. (I29j) . with the defini- 
tions 

c s = c+ , (30) 
we obtain the following recursion formula 

a2c's +3 + (2ai + 5a 2 )a 2 c' s+2 + (8a x + 10a 2 )c s+1 + (12ai + 10a 2 )c^' 
+alc s+3 + 2a 2 (a 2 + 2ai)c s+2 + (4(a^ + aia 2 ) + a 2 )c s+ i - 2a 2 c^ +3 
-(8a 2 + 6ai)c, +2 - (12a 2 + 18ai)c^ +1 - (8a 2 + ISaJc, = , (31) 
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and verify, for example, that the first two coefficients of the "+" series are 
related by 

+ _ (2a 2 (g 2 + 2^) + (2a! + 5a 2 )a 2 (7 + + 1)7+ - 7+(8a 2 + 6ai)) + 
ai + a 2 (7 + + l) 2 -a 2 ( 7+ + l) C ° 

For the "— " solution we just have to exchange c + for c~ and 7+ for — / y_ in 
Eqs.( l30|) to ( 132]) . If we keep only the leading terms in f±(x) we can write 
hf+lx) « 6ix" 7 +(l + 4/x) = and 6 a /_(ar) « 6 2 x 7 ~(l + c^/x) = 

firreg P - Applying the UV boundary condition we obtain that 



2 



V • ( 33 ) 



i(l + 7+ )+ 7+ )) 



The IR boundary conditions force the solutions to be equal to 1 as x — > 0, 
and the amount each solution (/ + or /_) contributes to the self-energy can 
be measured by the ratio R 



Masymp 

R oc Ueg 



rasymp 
J irreg 



(34) 

x~^A 2 /M 2 



The ratio R indicates which solution dominates the CSB. To gain some 
insight on the problem we can recall that when the gluons are massive the 
coupling constant (see Eq.((7])) is not expected to be larger than 0.5, therefore 
a 2 is a small number, and we can see that in the limit a2 ~ 0, the coefficients 
Ci and Ci are large and u = y/1 — 4a 2 m 1, leading us to 7_ = ~ and 



7+ = S*±l » 1, implying 



R 



2cJ 



(35) 

x^A 2 /M 2 



In the ultraviolet limit we verify that the ratio R behaves as (A 2 /M 2 ) u) and 
the f reg solution gives the asymptotic behavior of the self-energy. This is, 
apart from logarithmic contributions, the known 1/p 2 behavior found by 
Politzer using the operator product expansion [37j|. Notice that the asymp- 
totic behavior is fully described by the 1-gluon exchange, and this is not 



much different from what was found by Takeuchi [34j in a problem where 



the CSB is dominated by a four-fermion interaction. The influence of the 
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confining propagator enters only through the boundary conditions. As the 
mass solution that comes from the confining contribution has a fast falloff, 
almost behaving as an effective four-fermion interaction it is not surprising 
at all to see the similarity of our Eq. (l35j) to Eq.(29) of Ref. 34J. Another in- 



teresting limit of Eq.f l35j) is obtained when uj — >■ 0, or 4a2 = SC^g 2 /^ 71 " 2 ~ 1 ; 
although this may happen only for almost conformal theories, and it is a case 
of interest for technicolor models. We mention this because g 2 has an upper 
limit when the gauge bosons acquire a dynamically generated mass, as shown 
in Eq.Q j38|, and consequently we must force the j3 function coefficient b 
in Eq.([7|) to be very small in order to have co — > 0, and it may be quite 
difficult to build realistic technicolor theories with 4a2 ~ 1. This particular 
possibility is under study and will be presented elsewhere. 




4 0.00 



Figure 5: The ratio R = f^ mp '/ ' ft7^ P wh en A 2 /M 2 « m 2 /M 2 w 1 as a function of 
some natural values for the parameters a± and a2, which are proportional to the effective 
couplings for the different contributions entering into the gap equation. In this particular 
case the figure indicates that the irregular solution dominates over the regular one for a 
larger coupling value (oc ai) of the confining propagator. 



Equation ( |35l) would have a different (and interesting) behavior if the 
upper cutoff (A) were of order m ~ M. This would happen if the integration 
of the confining part of the gap equation, the one responsible for the CSB, 
were limited to a region in momentum smaller than, or of the order of the 
scale m. Assuming an upper limit in the momentum integration, we show 
in Fig.(jS]) the ratio R when A w m as a function of ai and a 2 . This ratio 
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is smaller than 1 indicating that the asymptotic behavior of the irregular 
symmetry breaking solution would dominate over the regular one in this 
particular limit. In order to explore even more this possibility we can assume 
that the confining contribution could be reduced to an effective four-fermion 
interaction. Some of the reasons why we are concerned with the possibility of 
generating a four-fermion interaction are the following: First, if confinement 
is introduced into the gap equation we should expect to reproduce some of 
the many phenomenological successful quark- mo dels based on the Nambu- 
Jona-Lasinio type of interaction. Secondly, lattice simulations show that 
the relevant gluonic energy scale of spontaneous CSB is due to the low- 
momentum component of the gluon field 3l|, |32|, |33j , which may indicate the 
possibility of a natural upper cutoff in the momentum. The existence of a 
specific momentum that separates the confinement and perturbative regions 
has also been discussed in a different context 39|. Finally, the existence 



of a completely nonperturbative infrared fixed point, as happens when the 
theory develops a dynamical gauge boson mass [25j], may induce effective 



four-fermion interactions as discussed many years ago in Ref . |40|, |41 . 

It is known that as long as we have a "massive" gluon propagator it 
would be possible to consider that this mass could be factorized from the 
propagator generating an effective four-fermion interaction, but this is not 
true because the actual interaction strength is measured by the product 
"coupling(g)propagator" , and we know from Eq. fl7|) that the 1-gluon exchange 
has not strength enough to generate such effective coupling. On the other 
hand the confining effective propagator, with the usual values for the string 
tension, is strong enough to generate the following effective gap equation: 



M 4f (p 2 



+ 



2 K F 
7r 3 m 4 

Co 



(2 



7T 



d 4 k 



d 4 k 



M 4f {k 2 



k 2 + M 2 f (k 2 ) 
9 2 (P 



6(m 2 - k 2 ) 
- k)3M 4f {k 2 ) 



[(p - k) 2 + m 2 g (p - k)][k 2 + M 2 f (k 2 )] ' 



(36) 



Apart from the gluon mass effect appearing in the 1-gluon contribution, the 



above equation has been extensively studied in Refs. 34j and 35[, and it 
does lead to a self-energy solution that decreases slowly with the momentum, 
or the so called irregular solution for the self-energy 35] . The solution of 
Eq. fl36l) follows from Refs. 3_4| and 35( observing the interplay between 
their 4-fermion coupling constant A and our effective coupling Kp/m 2 . 

The critical behavior of Eq. fl3"6"|) can be compared to the one of the com- 
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plete Eq. ffTBj) studied in the previous section, but now with a 4-fermion ef- 
fective kernel given by 



Kif\\ 2 = I dx I dy * I -^ J -^0(K-y) + 



oo px 



(I JO 

oo 



Jx 



y 2 I 


2K, 


(v + iy \ 




y 2 


{ *K f 







3C 2 g\x) 



16/T (x + e) 



We separated the kernel of Eq. (l37|l into two different kernels, one due to 
the effective four-fermion interaction and another due to the exchange of a 
massive gluon. Using the triangle inequality of Eq.( )20l) and the bifurcation 
condition given by Eq. (fl9l) we compute the critical condition and show in 
Fig-© the dot-dashed (black) curve of critical m values for the generation 
of massive solutions of Eq. fl36|) . This curve was obtained for m g = 600 MeV, 
rif = 3, Aqcd = 300 MeV and K F = 0.18 GeV 2 , and for comparison we also 
draw in Fig.()6]) the dashed (blue) critical curve of the complete kernel given 
by Eq. f[T7|) (without the four-fermion approximation) computed with the 
same parameters. This shows that most of the symmetry breaking is driven 
by the confining effective propagator and the four-fermion approximation is 
reasonable up to an order of 10%. Since we simplified the gap equation and 
used the triangle inequality to obtain the curve of Fig.©, it is difficult to say 
how much of the difference between these curves result from the numerical 
procedure. However it indicates that the confinement effect, as proposed in 



Ref . 19] , may indeed be the generator of an effective four-fermion interaction. 



We also computed the bifurcation condition (Eq. ffT9~|) ) with the part only 
due to the confining kernel in the four-fermion approximation for different 
values of the string tension. The result is shown in Fig. ([7]) where the con- 
tinuous (blue) curve was obtained for K F = 0.18 GeV 2 and the dot-dashed 
(black) and dashed (blue) were obtained respectively for Kp = 0.20 and 
Kp = 0.25 GeV 2 . Small changes in the string tension introduce larger effects 
than the ones that can be observed when we vary the gauge boson mass in 
the bifurcation case of the complete gap equation. This confirms again the 
dominance of the confining effective propagator. 

There is no doubt that the confining gap equation can be reduced to an 
effective four-fermion approximation. However we would like to claim that 
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Figure 6: The bifurcation condition given by Eq. (fT9|) in association with the triangle 
inequality of Eq. ([20]) is used to compute the critical m values for the massive solutions of 
Eq.©. The result for m g = 600 MeV, n f = 3, A QCD = 300 MeV and K F = 0.18 GeV 2 
is shown by the dot-dashed (black) curve. For comparison we also draw the dashed (blue) 
critical curve of the complete kernel given by Eq. (|17l) computed with the same parameters. 

the confining part of the fermionic SDE should have an upper cutoff at some 
scale not too much different from m. The reason for this is that the linear 
potential must break at some critical distance. For n/ = 2 quarks in the 
fundamental representation, lattice QCD data shows that the string breaks 
at the following critical distance [42| 

r c « 1.25 fm , (38) 

which corresponds to a m value compatible with the one necessary for the 
expected amount of CSB. 



The results obtained up to now contain many approximations, as the as- 
sumption that all mass scales were the same. A realistic calculation taking 
into account the different mass scales that we have in the original equation 
will be performed numerically in the next section, however we may see that 
the effect of the confining propagator dominates over the 1-gluon exchange 
and may even modify the asymptotic behavior of the gap equation as seen 
in the four-fermion approximation, which retains the CSB information. 



20 



0.05 0.10 0.15 0.20 0.25 0.30 0.35 
M(GeV) 



Figure 7: The bifurcation condition (Eq.([T9|)) is plotted only with the contribution of the 
confining kernel in the four-fcrmion approximation for different values of the string tension. 
The continuous (blue) curve was obtained for Kp — 0.18 GeV 2 and the dot-dashed (black) 
and dashed (blue) were obtained respectively for Kp = 0.20 and Kp = 0.25 GeV 2 . 



5. Chiral parameters with confinement and massive gluons 

In the above sections we studied aspects of the bifurcation equation as 
well as the asym pto tic behavior of the full gap equation, but to know how the 
proposal of Ref . [19| provides the right amount of chiral symmetry breaking 
it is necessary a full computation of the complete gap equation without the 



approximations made above or in Ref. 19j. We now solve Eq.flH]) without the 



angle approximation and taking into account the running coupling constant, 
with its infrared value dictated by a dynamical gluon mass, and a massive 
gluon propagator with the running mass included exactly as determined in 
the Ref. [6(. We choose the same values of n/, Aqcd, and Kp that were used 
to obtain Fig.©, but select an infrared value of the dynamical gluon mass 
equal to m g = 500 MeV, which is an average of many determinations of 
this parameter [5j. We compute numerically the full gap equation for the 
dynamical quark mass M(p 2 ) with different m values and show the results 
in Fig.©. 



A reasonable dynamical quark mass of (9(250) MeV is obtained with a m 
value equal to 180 MeV. This is a little smaller but totally consistent, within 
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Figure 8: Dynamical quark mass obtained with the numerical calculation of the full gap 
equation (Eq.®). 



the different approximations, to the one obtained in Ref. [19 . 

In order to show how reasonable the four-fermion approximation is to 
the complete gap equation we plot in Fig. © the numerical calculation of 
Eq.(l36l). This figure was obtained with m = 180 MeV, n/ = 3, Aqcd = 
300 MeV, Kp = 0.18 GeV 2 and m g = 500 MeV. The infrared value of 
the dynamical quark mass is overestimated by an O(40%) compared to the 
one of Fig.©, while it is also evident the slowly decreasing behavior of 
the dynamical fermion mass with the momentum as it was advanced in the 
previous section. The fact that the four-fermion approximation overestimates 
the dynamical quark mass was already observed in Fig.©. 



To confirm that the scenario of Ref. 19) is fully consistent with the 
CSB phenomenology, we compute several chiral parameters with the same 
values of rif, Aq C d, Kp and m g discussed above and m = 180 MeV, which 
leads to the usually expected dynamical quark mass of 250 MeV. The chiral 
parameters, computed in the Abelian gluon approximation, that we consider 
are: 

a) The pion decay constant /„., which is given at first order by the ex- 
pression |43|, [44 
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Figure 9: Dynamical quark mass obtained with the numerical calculation of the gap 
equation within the four-fermion approximation (Eq.(|36jl). 



f 2 

J TT 



dx 



xM(x) 



M(x) 



x dM(x) 
2 dx 



(39) 



4tt 2 J [x + M 2 (x)] 2 

There are many improvements to this formula, one of them determined in 
Ref.[45], which gives the following correction factor 



An 2 



doc oc 



x 



( dM(x) 
\ dx 



2[x + M 2 (x)f 



(40) 



and the pion decay constant /„- is obtained as the result of the sum 

fl = 11 + Sf* , (41) 



which is to be compared to the experimental value f w = 93 MeV 46 
b) The quark condensate (qq) at the scale /j 2 = 1 GeV 2 [l| 



M , xM(x) 
dx 



(42) 



4tt 2 J [x + M 2 (x)} ' 

which should be compared to typical value of the quark condensate {qq)(l 
GeV 2 ) = (229 ±9) MeV 3 0. 
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m 


/* 


u 


(qq) (lGeV) 


B 


[MeV] 


[MeV] 


[MeV] 


[MeV 3 ] 


[MeV 4 ] 


160 


62.0 


71.12 


169.03 


100.09 


180 


54.51 


61.83 


156.12 


84.96 


200 


46.59 


52.14 


142.0 


69.43 


Expected Values 


93 


93 


229 ±9 


146 



Table 1: Values of f ni (qq) and B obtained from Eqs. ([39|) -(|43 |) as a function of to. 



c) The MIT bag constant B [48 



B = 12tt 2 



In 



x + M 2 (x) 



x 



M 2 (x) 



x + M 2 (x) 



(43) 



which should be compared to the MIT bag constant value of 146 MeV 4 



49j, |50j, |51 



The results for the three parameters discussed above, computed as a func- 
tion of m, are shown in Table (pQ). Considering the simple rainbow approxi- 
mation for the gap equation we see that the results are below the expected 
values, but better choices for the vertex function as well as higher order 
corrections for these quantities will bring them closer to the experimental 
values. 



6. CSB for higher dimensional representations 

The study of CSB for fermions in higher dimensional representations is of 
interest because it is a possible way to verify how this mechanism is distinct 
from the confinement one, as well as it is important for technicolor model 
building. If a type of Casimir scaling as the one predicted by Eq.(jUj) occurs, 
we expect that for higher dimensional representations the CSB typical mass 
scale would be different from the one for the fundamental representation, and 
perhaps different from the confinement scale. It has also been argued that for 
"quarks" in the adjoint representation the dynamically massive gluons may 



have enough strength to generate a dynamical quark mass [181 . |20| . Indeed 
in Ref.[2o| a large dynamical mass was found for fermions in the adjoint 
representation, and we naively would expect that in this case the confining 
and chiral breaking transitions would appear separa tely . 



If we follow straightforwardly the model of Ref. [19( we must also verify 



what is the difference introduced by the confining propagator in the case of 
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higher dimensional representations, because in principle we should replace 
the string tension Kp by Kr, which is the string tension for fermions in the 
representation R. We assume that this replacement is accurate, although we 
know that the phenomenological potential of Eq.([T]), and consequently the 
string tension, does change according to the representation. For instance, in 
the case of the adjoint representation it is known that 

V A (r -> 00) = 2M 9 , (44) 

where M g is the energy of the lightest glueball. Moreover, the adjoint rep- 



resentation is not confined but screened [2JJ, what means that the confining 
propagator should be understood as effective up to a certain distance in 
these cases. Of course, no matter what the fermionic representation is we 
shall have a critical distance above which the string breaks. We assume that 
in the model of Ref . [19| most of the chiral symmetry breaking is still related 
to the form of Eq.(j3J), which does not get the chance to be probed at large 
distances, consequently we may still expect that most of the CSB is driven 
by the "confining" propagator. 

The fermion condensate is the most frequent quantity used to character- 
ize the chiral phase transition, and it is this quantity that we will analyze 
to investigate CSB for fermions in different representations. The fermion 
condensate is described by Eq. (l4"21) . It is easy to compute this quantity for 
different fermion representations if we consider the gap equation in the four- 
fermion approximation given by Eq. (l36p . perform the angle approximation 
and neglect the gauge boson mass in the propagator of the 1-gauge boson 
exchange contribution: 



,2 



o 



k 2 + M 2 f {k 2 ) 



+ (2tt)< J d p 2 [k 2 + M 2 f (k 2 )f [P k) 



(2tt) 4 J k 2 [k 2 + M 2 f {k 2 )} 

On the other hand we may write the fermion condensate for the representa- 
tion R in the following form 



= _ ^ r dx xM *w (46) 

11 A7r 2 J dX [x + M 2 R (x)] ' 
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where N R is the dimension of the fermion representation R and M R (x) its 
dynamical mass. We are forcing the upper limit of Eq. ( I46|) to be of 0(m), 
which can be as low as 0.18 GeV, while (qq) is well known at the 1 GeV scale. 
We did this, as we shall see below, in order to easily compare the condensate 
expression to Eq. fT4"5"]) . but we can check that Eq. (|4"B"]) provides a good estimate 
of the quark condensate at 1 GeV. To do so we see that in the four-fermion 
approximation the dynamical quark mass is almost flat up to 1 GeV, as 
can be seen in Fig.fjHJ), and its value is around 0.37, therefore the integral in 

Eq.flU]) can be approximated by /„* GeV2 cix[0.37x/(a; + 0.137)] « 0.263. With 
N R=3 = 3 for the fundamental representation we obtain (qq) = — ((0.27) 3 ) 
GeV 3 , which is a little large but consistent with the fact that the four-fermion 
approximation overestimates the dynamical mass. We stress that the upper 
limit in the first integral in the right-hand side of Eq. (145]) may be a physical 
one in order to be consistent with the critical distance at which the string 
breaks, as shown in Eq. (l38~l) . 

Eq.f l46p can be compared to Eq. fl45l) . the dynamical fermion mass M R in 
the four-fermion approximation, if we set all integrals at the scale m obtaining 



Since the first term between brackets in the right-hand side of Eq. fj47]) is 
much larger than the second one, combining the two last equations we have 



In the QCD case, for quarks in the fundamental representation, this relation 
underestimates the condensate due to the fact that the integration area in 
this equation was drastically reduced when we cutoff the integral at m 2 , 
whereas the dynamical quark mass solution, as shown in Fig.flH]), is almost 
flat up to 1 GeV 2 . Since the effect of the Eq.Q is the dominating one, it 
is quite plausible that the relation of Eq.flU]) holds up to other scales (still 
keeping the factor m 4 in the right-hand side) and it can be tested through 
lattice simulations. 

We can now make a few comments on the differences between CSB for 
fermions in the fundamental and adjoint representations estimating the ra- 
tio of the condensates for SU(3) fermions in these representations. First we 




(47) 



(qq) R (m 2 ) 



N R m 4 



M R (m 2 ) . 



(48) 
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need to know how the string tension changes as we change the fermion repre- 
sentation. It has been observed in lattice simulations what is usually called 
Casimir scaling for the string tension 21], i.e. 



K R « %K f , (49) 

Up 

where C^/Cp is the ratio between the Casimir operators for the represen- 
tation R and the fundamental one. For SU(N) theories and a finite N the 
Casimir scaling law must break down at some point, to be replaced by a 



dependence on the iV-ality k of the representation [21 



K R = f{k)K F (50) 

This change of behavior is credited to an effect of force screening by the 
gauge bosons. For fermions in the adjoint representation the iV-ality is zero, 
therefore, according to Casimir scaling, the adjoint string tension is given by 

2iV 2 

K A = j^Kf , (51) 

and, as a reasonable approximation, we may assume w IK p. Conse- 
quently we obtain the following ratio at the scale m 2 



(qq) 3 „ 3M 3 
(qq} 8 ~4M 8 



(52) 



Once the dynamical masses almost scale with the string tension value we 
could say that the above ratio is roughly of order 3/8. Of course, the un- 
certainty in this estimative is certainly connected to the remarks made at 
the beginning of this section about the phenomenological potential and the 
effective propagator for the adjoint representation. For other fermionic repre- 
sentations the screening behavior is smaller, although in all cases we certainly 
have a limit on the critical distance for which this approach is valid that will 
be connected with the string breaking mechanism. 



7. Conclusions 

In this work chiral symmetry breaking for QCD-like theories was studied 
in the case of a rainbow Schwinger- Dyson equation in the presence of dy- 
namically massive gauge bosons. Confinement was also introduced into the 
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SDE in the form of an effective propagator which is one that can reproduce 
an area law for quarks. This is a phenomenological way to investigate the 
possibility, indicated by the lattice, of confinement (by center vortices) being 
intrinsically related to chiral symmetry breaking. 

We briefly review the conditions for the confining propagator to generate 
non-trivial massive solutions. We studied the bifurcation condition for the 
complete gap equation, i.e. the one with the exchange of dynamically massive 
gauge bosons and with the inclusion of the confining effective propagator 
dependent on the entropic parameter m, which must be proportional to the 
dynamical quark mass, and on the string tension Kp, verifying that there is a 
maximum m value below which the chiral symmetry is broken, and generating 
the expected values for the dynamical quark masses. As already known in 
the literature we verified that the massive gauge boson exchange gives only a 
minute contribution to the dynamical fermion mass. Most of the breaking is 
due to the confining propagator and this may be one indication that we should 
not expect large differences between the chiral and confinement transitions. 

The asymptotic behavior of the gap equation was studied in order to eval- 
uate how confinement may affect the asymptotic self-energy. We observed 
that if the confining effect is restricted to a small momentum region, pre- 
senting some arguments for why this may happen, the solution changes its 
behavior from the one obtained with the help of the operator product ex- 
pansion to a self-energy typical of a four-fermion interaction. It is known 
that the massive gauge boson gap equation cannot be reduced to an effective 
four-fermion interaction due to the small infrared strength of the product 
"couplingCgipropagator" , however we verified that this is not the case for the 
confining effective propagator. The bifurcation equation for the gap equation 
with the four-fermion approximation performed in the confining part repro- 
duces approximately the result of the complete gap equation without any 
approximation. Actually this approximation overestimates the dynamically 
generated fermion mass by an 0(10%). 

The complete gap equation was computed numerically without any ap- 
proximation in Section V. Our result for the dynamical fermion mass is con- 



sistent with the one obtained by Cornwall [19j. We also computed several 
chiral parameters which resulted to be of the order of the expected exper- 
imental values. These parameters show dependence on the effective propa- 
gator scale m, however larger differences in the dynamical mass can also be 
observed with a small variation of the string tension value (see Fig.©). 
In Section VI we discuss the case of chiral symmetry breaking for fermionic 
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representations different than the fundamental one. We found a simple re- 
lation between the fermion condensate and the fermionic dynamical mass. 
This relation depends on the dimension of the fermion representation, m and 
Kp, and we assumed that the form of Eq.(j3]) still holds for different fermionic 
representations up to a certain distance. In principle such relation can be 
studied in lattice simulations. Since most of the symmetry breaking is due 
to the confining effective propagator and the string tension approximately 
scales with the Casimir operator, we do expect simple relations between the 
condensate values for different representations, but no matter what repre- 
sentation we choose it seems that in this scenario we shall not have large 
differences between the confining and chiral transition mass scales. 



can ac- 



Our results indicate that the CSB mechanism proposed in Ref. [19 
count for the expected values of several known chiral parameters. The model 
also seems to indicate that the CSB scale is connected to the confinement 
one even for fermions in higher dimensional representations. How far can 
we assume the four-fermion approximation discussed here for the purpose 
of practical calculations still needs further analysis. Finally, a more precise 
determination of the CSB observables can be obtained with the introduction 
of more sophisticated vertex functions and higher order corrections. 

Acknowledgment s 

We are indebted to A. C. Aguilar for discussions and help with the nu- 
merical calculation and are grateful to Prof. J. M. Cornwall for a clarify- 
ing exchange of correspondence and comments on a previous version of this 
manuscript. This research was partially supported by the Conselho Nacional 
de Desenvolvimento Cientifico e Tecnologico (CNPq) (AD and AAN), and 
Fundagao de Amparo a Pesquisa do Estado de Sao Paulo (FAPESP) (FAM). 

References 

[1] C. D. Roberts and A. G. Williams, Prog. Part. Nucl. Phys. 33, 477 
(1994). 

[2] A. Cucchieri and T. Mendes, PoS QCD-TNT 09, 026 (2009). 

[3] A. C. Aguilar, D. Binosi and J. Papavassiliou, Phys. Rev. D 78, 025010 
2008. 



29 



[4] A. C. Aguilar, D. Binosi and J. Papavassiliou, Phys. Rev. D 81, 125025 
(2010). 

[5] A. A. Natale, PoS QCD-TNT 09, 031 (2009). 

[6] J. M. Cornwall, Phys. Rev. D 26, 1453 (1982). 

[7] A. C. Aguilar and J. Papavassiliou, JHEP 0612, 012 (2006). 

[8] A. C. Aguilar and J. Papavassiliou, Eur. Phys. J. A 35, 189 (2008). 

[9] A. C. Aguilar and J. Papavassiliou, Phys. Rev. D 81, 034003 (2010). 

[10] D. Binosi and J. Papavassiliou, JHEP 0811, 063 (2008). 

[11] D. Binosi and J. Papavassiliou, Phys. Rept. 479, 1 (2009). 

[12] B. Haeri and M. B. Haeri, Phys. Rev. D 43, 3732 (1991). 

[13] A. A. Natale and P. S. Rodrigues da Silva, Phys. Lett. B 392, 444 (1997). 

[14] A. A. Natale and P. S. Rodrigues da Silva, Phys. Lett. B 390, 378 (1997). 

[15] P. O. Bowman, et al, Phys. Rev. D 78, 054509 (2008). 

[16] R. Hollwieser, M. Faber, J. Greensite, U. M. Heller and S. Olejnik, Phys. 
Rev. D 78, 054508 (2008). 

[17] P. O. Bowman, et al, Phys. Rev. D 84, 034501 (2011). 

[18] J. M. Cornwall, Invited talk at the conference "Approaches to Quan- 
tum Chromodynamics" , Oberwolz, Austria, September 2008, hep- 
ph/0812.0359. 

[19] J. M. Cornwall, Phys. Rev. D 83, 076001 (2011). 

[20] A. C. Aguilar and J. Papavassiliou, Phys. Rev. D 83, 014013 (2011). 

[21] J. Greensite, Prog. Part. Nucl. Phys. 51, 1 (2003). 

[22] P. Gonzalez, V. Mathieu and V. Vento, hep-ph/1108.2347. 

[23] W. J. Marciano, Phys. Rev. D 21, 2425 (1980). 



30 



[24] S. Raby, S. Dimopoulos and L. Susskind, Nucl. Phys. B 169, 373 (1980). 

[25] A. C. Aguilar, A. A. Natale and P. S. Rodrigues da Silva, Phys. Rev. 
Lett. 90, 152001 (2003); 

[26] A. C. Aguilar, A. Mihara and A. A. Natale, Phys. Rev. D 65, 054011 
(2002). 

[27] D. Atkinson, J. Math. Phys. 28, 2494 (1987). 

[28] D. Atkinson and P. W. Johnson, Phys. Rev. D 35, 1943 (1987). 

[29] D. Atkinson, V. P. Gusynin and P. Maris, Phys. Lett. B 303, 157 (1993). 

[30] C. D. Roberts and B. H. J. McKellar, Phys. Rev. D 41, 672 (1990). 

[31] A. Yamamoto and H. Suganuma , Phys. Rev. D 81, 014506 (2010). 

[32] A. Yamamoto and H. Suganuma, PoS Lattice 2010, 294 (2010), hep- 
lat/1008.1624; 

[33] H. Suganuma, et al., hep-lat/1103.4015. 

[34] T. Takeuchi, Phys. Rev. D 40, 2697 (1989). 

[35] K.-I. Kondo, S. Shuto and K. Yamawaki, Mod. Phys. Lett. A 6, 3385 
(1991). 

[36] F. W. J. Olver, Asymptotics and special functions, Academic Press, 
(1974). 

[37] H. D. Politzer, Nucl. Phys. B 117, 397 (1976). 

[38] J. M. Cornwall, Phys. Rev. D 80, 096001 (2009). 

[39] S. J. Brodsky and R. Shrock, Phys. Lett. B 666, 95 (2008). 

[40] W. A. Bardeen, C. N. Leung and S. T. Love, Phys. Rev. Lett. 56, 1230 
(1986). 

[41] W. A. Bardeen, C. N. Leung and S. T. Love, Nucl. Phys. B 273, 649 
(1986) 



31 



[42] G. S. Bali et al. [SESAM Collaboration], Phys. Rev. D 71, 114513 
(2005). 

[43] H. Pagels and S. Stokar, Phys. Rev. D 20, 2947 (1979). 
[44] J. M. Cornwall, Phys. Rev. D 22, 1452 (1980). 

[45] A. Barducci, R. Casalbuoni, M. Modugno, G. Pettini and R. Gatto, 
Phys. Lett. B 405, 173 (1997). 

[46] K. Nakamura et al. [Particle Data Group], J. Phys. G 37, 075021 (2010). 

[47] V. Gimenez, V.Lubicz, F. Mescia, V. Porretti and J. Reyes, Eur. Phys. 
J. C 41, 535 (2005) and references therein. 

[48] R. T. Cahill and C. D. Roberts, Phys. Rev. D 32, 2419 (1985). 

[49] A. Chodos, R. L. Jaffe, K. Johnson, C. B. Thorn and V. F. Weisskopf, 
Phys. Rev. D 9, 3471 (1974). 

[50] A. Chodos, R. L. Jaffe, K. Johnson and C. B. Thorn, Phys. Rev. D 10, 
2599 (1974). 

[51] T. De Grand, R. L. Jaffe, K. Johnson and J. Kiskis, Phys. Rev. D 12, 
2060 (1975). 



32 



